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I. INTRODUCTION 


The parton distributions in hadrons play a very important role in particle physics. Cross- 
sections for hard processes involving incoming hadrons can, based on the factorization theo¬ 


rems of QCD Q| , be calculated as a convolution of the parton distributions with parton level 
cross-sections calculated using perturbation theory. The parton distributions are universal 
in the sense that each hadron has a unique parton structure which can be used to calculate 
all hard processes involving that hadron. The normal interpretation of a parton distribu¬ 
tion fi{x,Q‘^) is as the probability to hnd a parton i (quark of some flavor or gluon) with 
a fraction x of the hadron momentum when probed by the momentum transfer Q^. The 
Q^-dependence is ve^ successfully described in perturbative QCD (PQCD) by the logQ^ 


3iy j 

i 


evolution equations which means that, given the input distributions in x at a scale Qq 
large enough for PQCD to be applicable, one can calculate the distributions at any higher 

However, this starting x-shape, which depends on non-perturbative QCD dynamics of 
the bound state hadron, has not yet been successfully derived from hrst principles. Instead 
they are obtained by fitting parameterizations to data, in particular structure function 


measurements in deep inelastic lepton-nucleon scattering (DIS), e.g. the CTEQ 
MRST ^ parameterizations. 


and 


Here we present further developments of our phenomenological model used to derive 
the parton distributions from simple assumptions regarding the non-perturbative properties 
of the hadron. The basic idea is to define the valence parton momentum distributions in 
the hadron rest frame where we assume that they are described by spherically symmetric 
Gaussians. The typical width of these distributions is a few hundred MeV from the Heisen¬ 
berg uncertainty relation applied to the hadron size. Sea partons are described through 
quantum fluctuations of the nucleon into baryon-meson pairs, having the same quantum 
numbers as the original nucleon, and the sea parton is identihed with a valence parton in 
the fluctuation meson. Through this mechanism, asymmetries in the non-perturbative sea 
parton distributions naturally arise, such as the difference between the u and d distributions 
and an asymmetry in the momentum distribution of s compared to s. 

The explanation of asymmetries in the nucleon sea by fluctuations of the nucleon into 


baryon-meson pairs also occurs in “meson cloud” models as reviewed 


in Q]. 


Although the 
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meson cloud models have an elaborate theoretical formalism (see e.g. j7|) one is forced to in¬ 
troduce phenomenological parts in order to obtain numerical results that can describe data. 
This concerns the choice of form factor in the nucleon to baryon-meson (iV — BM) splitting 
function and the parton distributions in both mesons and baryons. Our model invoke similar 
hadronic fluctuations as a basic quantum phenomenon, but differs from meson cloud models 
in important aspects. We use Gaussian momentum distributions both for hadrons in the 
fluctuations and for the partons in the hadrons, and the hadronic fluctuations are only con¬ 
sidered at the low starting scale Qq, where the starting parton densities are dehned and then 
evolved using standard perturbative QCD evolution. Thus, we have physically motivated 
parton momentum distributions in the nucleon as well as in the fluctuation hadrons, while 
meson cloud models usually use parameterizations. Our model results in a complete set 
of parton density functions, which we compare with different experimental data sets giving 
insights into different aspects of the non-perturbative nucleon structure. 

In section 2 we develop the details of the model and in section 3 we compare to the 
relevant experimental data as well as to conventional parameterizations of parton densities. 
Finally, section 4 gives a concluding discussion. 

II. THE MODEL 

A. Valence distributions 




Figure 1: Probing a valence parton in the proton and a sea parton in a hadronic fluctuation 
(letters are four-momenta). 


0 , givh 


This work is an extension of the previously presented physical model p, giving the 
momentum distributions of partons in the nucleon, as illustrated in Fig. ^ The model gives 
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the x-shape of the parton distributions of a hadron of mass M at a momentum transfer 
scale Qo ~ 1 GeV, i.e. xq{x,Ql) and xg{x,Ql), which provide an effective description of 
the non-perturbative dynamics of the bound state nucleon. Our approach is not intended 
to provide the full wave function for the hadron, but gives the four-momentum k of a single 
probed parton, whereas all other partons are treated collectively as a single remnant with 
four-momentum r (see Fig. P). In the nucleon rest frame, there is no preferred direction 
and hence the parton momentum distribution is spherically symmetric. The shape of the 
momentum distribution for a parton of type i and mass rrii is then taken as a Gaussian 


fi{k) = N{ai,mi) exp 


{ko-mi)^+kl+kl+kl \ 

2^1 ) 


( 1 ) 


where N is the normalization. The Gaussian form is not only a reasonable hrst guess, but 
can be motivated as the collective result of the many small momentum transfers affecting 
the parton through the non-perturbative bound state interactions which cannot be calcu¬ 
lated properly. The width cxj of the distribution should be of order hundred MeV from the 
Heisenberg uncertainty relation applied to the nucleon size, i.e. ai = l/d^. This Fermi 
motion inside the nucleon provides the ‘primordial transverse momentum’, which has been 
extracted from deep inelastic scattering data and found to be well described by a Gaus¬ 
sian distribution of a few hundred MeV width Q| giving phenomenological support for this 
description. 

The energy component does not have the same simple connection to the Heisenberg 
uncertainty relation. To keep the model simple and reduce the number of parameters, 
we assume a Gaussian distribution around the parton mass with the same width as the 
momentum fluctuations, such that partons can be off-shell at a soft scale of the binding 
interactions. This means a parton fluctuation life-time corresponding to the nucleon radius. 

The momentum fraction x of the parton is then dehned as the light-cone fraction x = 
k+/p+. Here, four-momenta are expressed as p = {p+,p-,p±) where the ‘plus’ and ‘minus’ 
components are p± = E pz and the ; 2 -axis dehned by the probe. The fraction x is 
then invariant under boosts along the ^-axis and equivalent to the conventional momentum 
fraction x = kzjpz in a frame where Pz is large (‘inhnite momentum’ frame). 

In order to obtain a kinematically allowed hnal state, one must impose the following 
constraints. The scattered parton must be on-shell or have a time-like virtuality (caus¬ 
ing hnal state QGD radiation), and an invariant mass less than the invariant mass of the 
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hadronic system. Furthermore, the hadron remnant r is obtained from energy-momentum 
conservation and must be time-like in order to project it on a hnal hadronic state through 
hadronization. Referring to Fig. ^ for the dehnitions of the momenta, one thus obtains: 

= {k + qY < = (p + qY (2a) 


r‘^ = {p — kY > 0 (2b) 

An important consequence of these constraints is to ensure that 0 < a; < 1 and f{x) ^ 0 as 

X —> 1. 

The parton distributions are obtained by integrating eq. ([Q) with these conditions. Using 
a Monte Carlo method this can be achieved numerically without approximations. For the 
simple case of the valence distributions it is possible to derive an analytical expression as 
well, see sec. IIIBI below. The normalization of the valence distributions is provided by the 
sum rules 

/ dx Uy{x) = 2 and / dx dy{x) = 1 (3) 

Jo Jo 

to get the correct quantum numbers of the proton (and similarly for other hadrons). The 
normalization of the gluon distribution is chosen so as to saturate the momentum sum rule 

V /" dxxfi{x) = l. (4) 

where the sum is over the valence quarks and gluons. 


B. Analytical expressions 

By integrating f{k) in eq. over k\ and k- of the parton between the limits given by 
m, we arrive at the following analytical expression for the bare hadron valence distributions: 


fi{x) = N{ai,Xi) 
1 + erf 


1 -|- erf 

2di 


I- Xj 
2di 


exp 


(X - Xi 


4d? 


exp 


(i-x,Y 


( 5 ) 


where cTj = Xt = rrii/M and erf(x) = ^ Jq (fh® error function). The normal¬ 

ization N{di,Xi) is given by Q- This analytical expression can be simplihed when cxj and 
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rrii are sufficiently small {e.g. ai, rrii < M/4) to 


fi{x) 


N{ai,Xi) 


1 — exp 


f (1 

V 2a/ 



{x - a;i)2\ 

4a/ ) 


( 6 ) 



Figure 2: Valence u-quark distribution of the proton obtained from the model: the “bare” proton 
distribution xu^q{x) from Monte Carlo simulation compared to the analytical form eq. 0, as well 
as the “full” Ut,-distribution including simulation of partons in baryons in hadronic fluctuations 
Ip) —> \BM) as described in sec. Ill Cl below. The dotted line is the valence distribution of Q], 
having a lower width cr^ = 180 MeV but no effect of baryons in fluctuations. (All curves show the 
momentum density, whereas the normalization is in the number density according to eq. ©.) 

The “bare” proton valence u-quark momentum distribution xuQy{x) given by eq. © is 
shown by the full curve in Fig. |21 The approximation in eq. (ED gives an indistinguishable 
result and is therefore not shown explicitly. Fig. |21 also demonstrates a very close agreement 
with the Monte Carlo simulation result of the model (histogram). 

This behavior of the valence distributions is modihed by the contributions from the 
baryons in the hadronic fluctuations described below. At a: > 0.7 the shape of the valence 
distributions are still given exclusively by the original hadron, i.e. by eq. but at lower 
X the valence distributions are significantly modified. The shift of the distribution to lower 
X is clearly seen in Fig. |21 as the difference between the “bare” xuov{x) distribution and 
the “full” xUy{x) valence u-quark momentum distributions. A similar end result can also 
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be obtained without including this baryon fluctuation effect, as in P|, by a lower Gaussian 
width when fitting the model to F 2 data (discussed below in section IIII Afl . 

One can note that the asymptotic behavior for very large x > 0.9 is (1 — according 
to eq. ©• For 0.3 < a; < 0.9 however, the u quark distribution is well described by a form 
u{x) = 8(1 — a:^'^, in good agreement with the (1 — a:)^-dependence predicted by spectator 


counting rules 


Q. 


C. Hadronic fluctuations and sea distributions 

Sea partons arise from the non-perturbative dynamics of the bound state hadron, for 
which it should be appropriate to use a hadronic quantum mechanical basis. Therefore we 
consider hadronic fluctuations of the proton 

\p) = «obo) + OpTrlpTF^) a„^|n7r+) ... -f aAxl-^K^) + ... (7) 

where the different states are assumed to be orthonormal and the normalization constants 
are real and fulfill Og + J2bm ^bm = 1- 

Probing a parton i in a hadron H of such a fluctuation (Fig. gives a sea parton with 
light-cone fraction x = xnXp of the target proton, i.e. the sea distributions are obtained 
from a convolution of the momentum K of the hadron and the momentum k of the parton in 
that hadron. We assume that the momentum K of the probed hadron is given by a similar 
Gaussian as in eq. o, with a separate width parameter an (for simplicity and definiteness 
we make this parameter common for all fluctuations). The momentum K' of the other 
hadron in the fluctuation is then fixed by momentum conservation in the rest frame of the 
original nucleon. We let both hadrons be on-shell, which Axes their energies. This implies 
that energy is not conserved at this intermediate stage, but is of course restored for the 
observable final state. With the hadron four-vectors specified, one obtains the light-cone 
fraction xh = 7f+/(iF -|- iF^)+. A few possible modifications of the details of the model are 
discussed in the next section and the possible use of an effective meson mass instead of the 
physical pion mass is discussed in sec. IIII Gl 

The above model for valence distributions is then applied to the fluctuation hadron H to 
get the parton momentum and light-cone fraction Xp = in H. The flavor sum rules 

in eq. © must, of course, be modified to apply for H. The kinematical constraints of eq. 
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m are in this case modified to 


mi<f<W^^ = {K + qY (8a) 

> 0 and (r + = {p — > 0 (8b) 

such that the scattered parton and the remnants (c/. Fig. ^) are on-shell or have positive 
virtualities within physically allowed limits and hence can give a proper final hadronic state. 

A Monte Carlo method is used to simulate this two-step process by choosing K and fc, 
impose the constraints and obtain the momentum fraction x. By iterating the procedure the 
additional quark and gluon distributions due to the fluctuations are generated. Note that 
the flavor number density dx fi{x) is not affected by the convolution with the fluctuation 
momentum. This means that the normalization criteria 0 (modified for the hadron in 
question) automatically ensure flavor conservation, since the meson and the baryon in the 
fluctuation are multiplied by the same normalization constant 

Since the x distribution of the fluctuation hadrons is given by xh = K^/{K -|- K')+ ^ 
Mh/{Mb + Mm), the low-a; sea is mainly given by scattering off the meson in the baryon- 
meson fluctuation. The scattering off the baryon in the baryon-meson fluctuation gives 
a contribution resembling the valence distribution, but at lower x than the bare hadron 
distribution. This means that the analytical valence distribution of the bare hadron derived 
above is modified for x < 0.7 with the inclusion of hadronic fluctuations, as shown in Fig. El 

The normalization of the sea distributions is given by the amplitude coefficients cx\j^. 
These partly depend on Clebsch-Gordan coefficients and should also include a suppression 
for larger masses. However, the full dependence is due to non-perturbative dynamics that 
cannot be calculated from first principles in QCD, including the possibility of mixing of 
different mesons/baryons with the same quark content. We therefore take the normalizations 
^‘bm parameters. 

1. Alternatives in the definitions of fluctuation momenta 

In the details of the model described above we have made some choices which could be 


done differently: 


1. In the definition of the fluctuation momentum, the fluctuation hadrons were made on- 
shell. Another possibility would be to let the energy of the probed hadron fluctuate 
around the hadron mass as in eq. (^Q), allowing the fluctuation hadrons to be off-shell 
on a soft scale. The differences in the resulting sea distributions are however small 
and can be mimicked by changing parameter values. 

2. We defined the upper bound on the scattered parton virtuality as < Wjj = {K + 
qY- This choice keeps the parton distributions in the fluctuation hadrons close to 
unaltered with respect to the distributions in the corresponding free hadrons. However, 
this choice allows the struck parton to have essentially all the proton momentum, 
which might seem unnatural. An alternative approach would be to put the upper 
bound on to {xhP + qY — xhW^, such that it is restricted by the longitudinal 
momentum fraction taken by the fluctuation hadron. This results in a softer meson 
momentum spectrum, because the distributions in the fluctuation meson is distorted. 
For completeness, we will consider both alternatives in our predictions for the s — s 
asymmetry (sec. ITTTdI) . 

3. The limits on the remnant momenta was put to > 0 (to ensure that f{x) —> 0 as 
X —>• 1) and {r + RY = {p — kY > 0. This latter choice is motivated by the fact that 
the fluctuation hadrons are close in space, such that an exchange of momentum on 
a soft scale should be possible and therefore the condition is on the whole remnant 
system. Another possibility would be to force r and R separately to be time-like. This 
is obviously not possible if we scatter on a baryon with mass larger than the original 
hadron mass and make the fluctuation hadrons on-shell, since R^ > 0 can then never be 
satisfied. In the off-shell case discussed under point 1 this is possible, and gives slightly 
harder meson spectra for low meson masses, whereas for meson masses > 300 MeV 
the difference is negligible. 

D. GVDM contributions at low 

In the same spirit as the hadronic fluctuations of the proton described above in eq. dZD, 
one should also note that a photon may appear as a vector meson, such that the quantum 


9 
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Figure 3: Schematic view of the photon fluctuating to a vector meson before interacting with the 
proton. 


state should be expressed as 

I7) = C'olbo) + 


fv 


■\v) 


dm{- ■ ■ 


(9) 


y j y Jmo 

In the original vector meson dominance model (VDM) only fluctuations to vector mesons 
were considered j^, but in generalized models (GVDM), also a contribution from a con¬ 
tinuum of higher-mass states is included represented in eq. 0 by an integral over 
masses. 

The GVDM interpretation of the photon-proton interaction is illustrated in £g. |21 The 
corresponding cross section is given by a convolution of the photon-to-meson fluctuation 


probability with the meson propagator and the meson-proton cross section 


lll |. and the 


continuum contribution is included using a phenomenologically chosen spectral weight func¬ 


tion. As shown in 


F2{x,Q^) = 


12j |. the resulting expression for the proton structure function, 
(1 — x)Q‘^ \ f ^ ^ 


dvr^Q! 




rv 


- V=p,oJ,<f> 


V 


Q2 + 


1 + 


m 


V 


+ rc 


(l-ec); 


run 


+ ec^ln(l + ^) 




m. 


0 J 


Q 


2e 


A-,^ . (10) 


X' 


Q2 + 77^2 

describes data very well for < 0.7 GeV^. 

At higher this GVDM contribution should be phased out in order to conform to the 
conventional description in terms of parton density functions. One way of doing this Q 
is by introducing a phenomenological form factor I'he GVDM contribution 

for which gives a good description of HERA F 2 data at these intermediate 

< 4 GeV^. The GVDM component was found to be negligible for > 4 GeV^. 

The details of the phasing out of the GVDM component is not known and would require 
additional assumptions and parameters in the model. Fortunately, GVDM (giving symmetric 
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quark-antiquark contributions) can be neglected in this study of the asymmetries in the 
proton, since the relevant data used in the following is at > 4 GeV^ where the GVDM 
component was found to be negligibly small. 


III. RESULTS AND COMPARISONS TO EXPERIMENT 


With this simple model, based only on Gaussian momentum fluctuations and hadronic 
quantum fluctuations together with kinematical constraints, we are able to describe several 
different aspects of the parton distribution functions of the nucleon. Perhaps phenomeno¬ 
logically most interesting, we describe the asymmetries between u and d in the proton quark 
sea. Also, we get an asymmetry in the nucleon strange sea, such that the s quark has a 
harder momentum distribution than the s antiquark, which is of interest Q in relation to 
the NuTeV anomaly [l^ . 

Since we have most data on the proton, the parton widths are htted for this special 
case of a hadron. For other hadrons, we assume for simplicity that the widths of the gluon 
distribution is the same as that for the proton. For valence quarks, we assume that the d 
quark width is the typical width for a quark with quark number J q{x)dx = 1 in a hadron, 
while the u quark width is typical for a quark with number J q{x)dx = 2. This means 
that in all mesons and also e.g. the A baryon, all valence quark widths are given by the 
proton d quark width. This can be motivated by the hadrons having essentially the same 
spacial extent, making the only relevant difference the possible effects of the Pauli exclusion 
principle. Such ’Pauli blocking’ (see and references therein) would reduce the effective 
space available for the u quarks in the proton, resulting in a larger momentum fluctuation 
width. Therefore, for the neutron we simply make an isospin transformation d ^ u oi the 
proton distribution parameters. 

The model provides valence and sea parton x-distributions as shown in Fig. ^ These 
apply at the low starting scale Qg, while the parton distributions and the proton structure 
function F 2 {x, Q"^) at higher are obtained by applying standard QGD evolution as imple¬ 
mented in the QGDNUM16 code For the results below we have used the conventional 


MS scheme and the evolution equations in next-to-leading order to have the highest avail¬ 
able theoretical accuracy. Using the option of only leading order equations, the values of 
all model parameters are unaffected except Qq and ag as discussed in section 1111 El With 


11 







(a) (b) 


Figure 4: Parton distributions from our model at the starting scale Q^: (a) all partons distributions 
(b) the sea parton distributions (note the different scales). 


leading order evolution, however, the quality of the £t to data becomes less good. 

Since different model parameters are sensitive to different data, we have used data from 
several different experiments: 


Fixed-target F 2 data [lallll to £x large-a: (valence) distributions 

HERA F 2 data [l^ for the gluon distribution width and the starting scale Ql 


charge asymmetry data jl9| as a cross-check on the ratio of Gaussian widths for 
the u and d valence quark distributions 

J/h-asymmetry data 0 for the normalizations of the \pn^) and \n7T~^) fluctuations 


Strange sea data |21! to £x the normalization of the fluctuations including strange 
quarks. 


It should be noted, however, that most of the model parameters influence several observables 
at least to some degree, making the task of a total fit quite involved. 

It is interesting to note that this simple model can describe such a wealth of different 
data with just one or two parameters per data set, as will be seen in detail below. The 


12 













X X X X 


Figure 5: 


The proton structure function F 2 {x,Q'^); HI data 



compared to our model with its 


sensitivity to a ±50% variation of the width parameter ag of the gluon distribution. 


best-fit parameters are: 


au = 230 MeV ad = 170 MeV ag = 77 MeV 
Qo = 0.75 GeV an = 100 MeV (11) 

< = 0.45 < = 0.14 = 0.05 


The normalization for the proton fluctnating to a baryon B and a meson M is 

defined as the fraction of the proton momentnm taken by all partons i in the fluctuation 


BM, 



^xfi/six) + ^xfi/M{x) 
. ieB i&M 


( 12 ) 


A. Inclusive ±2 data 

The model was originally P] used only to fit inclusive data, using a somewhat different 
definition for the hadronic fluctuations, and was shown to give a good fit. Now we are mainly 
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figure 6: The proton structure function F 2 {x,Q‘^) for large x values; NMC and BCDMS data 
compared to our model, also showing the results of ±20% variations of the width parameters 
au and ad for the u and d valence distributions. 


interested in the model’s description of asymmetries in the proton, and hence the F 2 data 
is used mainly to £x some model parameters which are not well restricted by asymmetry 
data, namely the Gaussian width of the gluon distribution, cr^, and the starting scale Qq. 
However, in view of the simplicity of the model, we still get a nice description of F 2 data 
from HI j^, see Fig. El To show the sensitivity to the value of ag, we also show the result 
of varying this parameter with ±50%. The surprisingly low value for ag obtained in the £t 
will be discussed in sec. null 

Since this data is restricted to relatively low x values, we also compare to fixed-target 
NMC and BCDMS data HQ, wMch .. Ga..an wU.ha of the a aad d ,aaoh 

distributions, au and ad- This is shown in Fig. El which also shows the sensitivity of the £t 
to changing these parameters with ±20%. As seen in the figure, the u width is quite well 
constrained by the large-x data, while the d width is less well constrained, and the data 
might account for a ad which is up to 20% smaller. This is certainly true for the NMC 
data, which actually seems to prefer a smaller d width. However, the full curve describes 
the best-fit value due to the very small error bars on the BCDMS data, so we keep this as 
our main result in the following. 
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Figure 7: The charge aswnmetry versus rapidity, eq. d, for leptons from hF^-decays in pp 
collisions at the Tevatron 1^ (with positive and negative lepton rapidity combined) compared to 
PYTHIA simulations using parton densities from our model, with best-fit parameters and a 20% 
reduced width of the valence d quark distribution. 


B. The forward-backward asymmetry 


The forward-backward charged lepton asymmetry from decays in pp collisions at 
the Tevatron provides information on the large-a: distributions of d and u quarks. This is 
because the W bosons are produced mainly through the processes ud —> W~^ and du —>• W~. 
If the u distribution is harder than the d distribution, we will get more W~^ than W~ in 
the direction of the proton beam, and vice versa in the direction of the p-beam, which is 
precisely what is observed. The data is for and 0.006 < x < 0.34, corresponding 

via PQCD evolution to 0.01 < a: < 0.6 at the starting scale Qo ~ 1 GeV^. We can therefore 
use the data on the charge asymmetry versus rapidity, expressed as 

da^/dyi — da~ / dyi 


Mvi) = 


da~^/dyi -|- da /dyi 


(13) 


to get further information on the Gaussian width of the d quark distribution. The result, 
using our PQGD-evolved parton distributions in PYTHIA to simulate production via 
the processes qq' —>• W^, qq' —>• gW^ and qg —*■ W^q', and subsequent decay 
is shown in Fig. [Tj It is clear that the model reproduces the salient features of the data, in 
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particular when the d width is reduced as preferred by the NMC data discussed above. 

The difference between the widths of the u and d quarks ((j„ = 230 MeV and ad = 135 — 
170 MeV) can perhaps be attributed to the Pauli blocking effect on the u quarks as discussed 
above. Large differences in the intrinsic transverse momentum between the u and d quarks 
are indeed indicated by data from HERMES j^. Also the standard parameterizations of 


parton densities 


result in harder spectra for the u than for the d valence distribution. 


C. The d - u asymmetry 

Now we come to the important results of the model regarding the sea quark asymme¬ 
tries. These result from the fluctuations of the proton to meson-baryon pairs. The hrst 
such asymmetry that we will investigate, and the one with most experimental data, is the 
difference between u and d distributions. 

In this model, the u and d sea comes mainly from fluctuations of the proton into a nucleon 
and a pion, since these are the lowest energy fluctuations. This means that we will get an 
excess of d over h, since |pvr°) is symmetric between d and h, while |n7r+) contains only d and 
no u. If these were the only available fluctuations we would get from isospin Clebsch-Gordan 
coefficients the ratio = 1/2. However, this is modihed by the inclusion 

of IA7r)-fluctuations (which have small mass suppression relative to the |V7r) fluctuations). 
In the SU(6) scheme, the IA^+tt”) fluctuation has a larger Clebsch-Gordan coefficient than 
the IA'^tt’*'), which increases the relative amount of u quarks. In our model, we take such 
higher mass fluctuations into account by simply htting the normalization parameters 
and using experimental data on the light sea. 

The ratio d/u can be extracted by comparing muon pair production through the Drell-Yan 
process, qq , in pp collisions with that in pd collisions assuming isospin symmetry. 

Using parameterizations for the sum d{x) + u{x), also the difference d{x) — u[x) can be 
extracted. 

We have used data from the E866 collaboration Q to £t the parameters relevant for 
the d and u distributions. These parameters are the normalizations and and the 
Gaussian width an- However, the sensitivity of the meson parton distributions to an is 
small, and in particular the position of the peak of the momentum distributions is quite 
independent of an, and depends only on the mass of the fluctuation hadrons. For small 
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Figure 8: Comparison between our model and data from the E866/NuSea collaboration 


20|: (a) 


u{x)/d{x) (b) d{x) — u{x). The full line shows the results of the model using the physical pion 
mass, while the dashed line shows the results with an effective pious mass = 400 MeV as 
discussed in the text. The parton distributions are evaluated at = 54 GeV^. 


an ^ 100 MeV, also the overall shape of the meson distributions is mainly determined by 
the hadron masses. For larger an ^ 100 MeV, the parton distributions inside the fluctuation 
hadrons are altered due to the kinematical limits in such a way as to keep the convoluted 
shape stable. Therefore we take the value of an to be 100 MeV, which is preferred by the 
fit. 

In the distribution of d{x) —u{x) (or, equivalently xd{x) —xu{x)), the effect of symmetric 
qq pair production from gluons in the PQCD evolution is canceled and thus the result 
depends only on the properties of the non-perturbative distributions at the starting scale, 
barring the uncertainties from the extraction of the data from experiments. Therefore it is 
quite interesting to compare the experimental results for this distribution with the model 
results. As can be readily seen in Fig.|H|D, the pion distribution using the physical pion mass 
tUtt = 140 MeV (full line) is too soft, peaking at x ~ 0.05, while the data peaks around 
X ~ 0.1. This cannot be cured using a larger cth , as explained above, but rather indicates 
that one should use a larger effective pion mass. Allowing this effective pion mass to be a 
free parameter in the fit, the best-fit result is that rrf^ ~ 400 MeV, as showed in Fig. |H1 
(dashed line). Note that the points driving the £t are the small-x d/u data, which have the 
smallest errors. 
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One can think of several explanations for such a large effective pion mass. One could 
be that we need to include a relatively large fraction of heavier mesons, such as the p 
mesons, in our description. This would be in accordance with the finding of a large \Np) 
component in the study based on a meson cloud model. Another explanation would be 
that since the pions have very non-typical hadron masses, the meson-baryon fluctuations 
involve more generic meson states, with the quantum numbers defined by the actual hadron 
states available, but with masses with some distribution which can be approximated with 
a meson mass of 400 MeV. In any case, the introduction of a larger pion mass without 
any other modifications of the model reduces the quality of the fit to the low-x F 2 data 
described in sec. 1111 Al (the effect is comparable to changing the value of ag by less than 
-|-50%), which points to the need of introducing more assumptions and parameters to get 
a perfect description of all data. Since our intention is not to get the best possible fit (as 
conventional parton density parameterizations), we refrain from doing this without a physics 
motivation related to our model. 

D. The strange sea asymmetry 


Another consequence of the model is an asymmetry between the strange and anti-strange 
quark ^-distributions in the nucleon, an aspect that we did a special study of in [l^ and 
extend here. The reason for this asymmetry is that in fluctuations with open strangeness, 
the s quark is found in the baryon while the s quark is in the meson. Since the baryon, due 
to its higher mass, has a harder a;//-spectrum than the meson, we get a harder x-spectrum 
for the s quark then the s. This effect persists in spite of the fact that the s distribution in 
the A peaks at lower Xp than the s distribution of the see P| for examples of typical 
meson and baryon distributions from the model. 

For simplicity we assume that all fluctuations including strangeness (such as \ AK*), |SiF)) 
can be implicitly included in the |AA')-fluctuation. Fluctuations where the ss pair is part 
of a meson wave function, thus giving a symmetric contribution, are neglected since they 
should be suppressed (0 due to large mass and p due to Clebsch-Gordan coefficients). The 
resulting s and s distributions are shown in Fig. IHd. 

Note that in this model the criterion that the nucleon should have zero total strangeness, 
Jq (s(x) — s(x)) = 0, is automatically fulfilled, since all strange fluctuations have the same 
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Figure 9: CCFR deep inelastic scattering data j21(] on the strange sea distribution (rEs(x) + 
xs{x))/2 in the nucleon at different compared to our model based on \A.K) ffuctuations (with 
normalization ‘s/u’= (s + s)/('U + d) as discussed in the text) at a low scale and evolved to larger 
with perturbative QCD that also adds a symmetric perturbative ss component. 


number of s quarks in the baryon as s quarks in the meson due to the sum rules Q modified 
for the hadrons in question. This normalization is not changed when the parton distributions 
of the fluctuation are convoluted with the distributions of the fluctuation hadrons in the 
nucleon (except for the factor asM which is common to the baryon and the meson). 

Since the normalization of the \AK+) fluctuation cannot be safely calculated, we take 
it as a free parameter which we £t to data from the CCFR collaboration 0 on the av¬ 
eraged strange sea (s -|- s)/2, as shown in Fig. IHl Here, perturbative QCD evolution to 
larger shifts the original s and s distributions to smaller x and adds a symmetric 
sea arising from g —^ ss. The fit gives a |AiF-|-) normalization = 0.055 such that 
Jq dx{xs{x) + xs{x))/ dx{xu{x) -f xd{x)) ~ 0.5, i.e. the strange sea momentum fraction 
at Qq is approximately half of that of a light sea quark, in agreement with the parton density 
analyzes in e.g. 0 and the meson cloud model study in 

This normalization of the non-perturbative strange quark sea means that the coefficient 
has essentially been scaled down with the fluctuation time At ~ 1/AE relative to 
the light sea. This is a smaller suppression than 1/(AF^)^ that one might have expected 
based on old-fashioned perturbation theory. On the other hand, this ss suppression is of 
similar magnitude as the ratio P{ss)/P{uu) ~ 1/3 of the probabilities for quark-antiquark 


production in phenomenological hadronization models, such as the Lund model 


2J. Given 
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X 


Figure 10: The strange sea asymmetry s“(x) = xs{x) — xs{x) (at = 20 GeV^) from the 
model and combined with the function F{x) accounting for NuTeV’s analysis giving 6w = 

Jq dx s~ {x)F{x) = —0.0017. 

that both cases concern non-perturbative ss pair production in a color held, this need not 
be surprising but indicate common features. 

This \AK+) huctuation results in an asymmetry S~ = dx{xs{x) — xs{x)) = 0.00165 
at = 20 GeV^. However, this result is not altogether independent of the details of the 
model. The precise dehnition for the upper limit of the struck quark virtuality discussed in 
sec. HTTTTl item 121 and the width of the d distribution (which gives the width of the valence 
quark distributions in the A and K) both have inhuence on the value. Using both dehnitions 
for the limit, and allowing the d width to vary between the values used in Fig. [3 we get a 
range of values for the asymmetry 


0.0010 < ^“ < 0.0023 


(14) 


An asymmetry in the momentum distributions of s and s is especially interesting in 
view of the NuTeV anomaly. In the NuTeV experiment the value of siv? 9w was 
extracted from neutral and charged current cross-sections of neutrinos and anti-neutrinos. 
The value they hnd differs by about 3cr from the value obtained in Standard Model hts to 
data from other experiments: sin^ = 0.2277±0.0016 while sin^ 6^ = 0.2227±0.0004. 


A number of possible explanations 


25j | for this discrepancy have been suggested, both in 


terms of extensions to the Standard Model and in terms of effects within the Standard Model. 
One explanation within the Standard Model is that, since the neutrinos and anti-neutrinos 
interact differently with s and s quarks, an asymmetry S'" 7 ^ 0 shifts the NuTeV result. In 
order to facilitate the calculation of such a shift, NuTeV has published a folding function 
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F{x) to account for their analysis and give the shift in the extracted value of siv?9w^ 
where the parton densities are to be calculated at = 20 GeV^. Using this, we get the shift 
Asin^6*vi/ = dx s~ {x)F{x) = —0.0017, corresponding to a reduction of the significance 
of the anomaly with one standard deviation. Our results for s~{x) = xs{x) — xs{x) and 
s~{x)F{x) are shown in Fig. The range of values of S~ of eq. (d corresponds to 
—0.0024 < Asin^^w < —0.00097 for the shift in sin^ This reduces the anomaly to 

between 1.6a and 2.4cr away from the Standard Model value. 

This important topic has been investigated previously. On the theoretical side we note 
that it has recently been shown 22| that higher order perturbative effects give a negative 
contribution to S~, although significantly smaller than our positive non-perturbative effect. 


Within the meson cloud model of 


28j |. based on exponentially suppressed AK fluctuations. 


an inconclusive result on the sign of S was obtained and a strange sea that was too small 


to affect the NuTeV anomaly. The meson cloud model 


in 0] 


obtained that the s quark 


has the softer spectrum than the s, corresponding to a negative S~. Using a light-cone 
two-body wave function model applied to AK fluctuations |2^, a positive result similar to 
ours was obtained in H. and . 0 an even lan.e. eednetion of fde nno.nnl, fKnn on. wne 
found using an effective chiral quark model. Unfortunately, none of these studies provide 
comparisons with the measured strange sea, making the signihcance of their results difficult 
to assess. 


The experimental situation is at present unclear. 


Ln 


. 0 . 


a positive S was favored based 


on 


several earlier experiments, whereas NuTeV Lm obtains the opposite sign based 


on 


their own data. However, the global analysis in of the s — s asymmetry, including the 
NuTeV data as well as the CCFR data and using a very general functional form for s“(x), 
gives a best fit value for the asymmetry S~ of the same magnitude and sign as ours. 

It is clear that an asymmetry in the strange quark sea of the nucleon arises naturally in 
models where non-perturbative sea quark distributions originate from hadronic fluctuations 
of the nucleon. Our model, which reproduces the measured strange sea, reduces the NuTeV 
anomaly to a level which does not give a significant indication of physics beyond the Standard 
Model. 
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E. Comparison to standard parameterizations 


Let us first note that it is a non-trivial result that the many different data sets shown 
above can be well described by a model with only the eight parameters in eq. (Hill . The 
model provides all the different parton distributions in x at Qq, he. 9, u, d, 

s, s. This is achieved with only four shape parameters (the cr’s) and three normalization 
parameters (the he. signihcantly fewer shape and normalization parameters than 

distributions. 




The simplest possible conventional parameterization based on counting rules M: xfi{x) = 
Aix{l — x)“% cannot be as economic in the number of parameters and still allow different 
shapes of different sea distributions. Considering only valence and light sea distributions 
and neglecting the strange sea, we were able to get reasonable hts of the same data sets 
with eight free parameters (when applying the constraints in eq. ()dl4|l i. However, the htted 
values for the exponents Oj for the sea quarks and the gluon turn out more than twice as 
large as expected from the counting rules (using the expected values gives a poor description 
of the data). With our model we obtain a generally better description of the data, which is 
achieved with fewer parameters. Moreover, a physically motivated model whose parameters 
have a physical interpretation can give new insights about the unknown non-perturbative 
QCD dynamics embodied in the parton distributions. 


Q. 


The standard parameterizations of parton densities, such as CTEQ m and MRST 
have many more parameters (typically 20 for the x-shapes and additional ones for normal¬ 
izations) and provide high-quality hts to large sets of different kinds of data. It is, therefore, 
interesting to compare the parton densities resulting from our model with such standard 
parameterizations. 

In Fig. [TT] we show our parton densities together with the uncertainty band of the 


CTEQ6M distributions 
are striking: 


at the low scale = 1.3 GeV^. Several features of the plot 


1. For large and intermediate x, our valence quark as well as sea quark distributions from 
the model agree quite well with those from CTEQ6M. 


2. Our gluon distribution is slightly lower at large x and much larger at small x. 

3. For small x < 10“^ the quark and anti-quark distributions of CTEQ are clearly higher. 
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(b) 


Figure 11: Comparison between parton densities at = 1.3 GeV^ from our model (thick lines) 


and the CTEQ6M uncertainty band (plotted using the Durham HEPDATA server 


^1). Note 


the different scales on the axes of (a) and (b) and that the gluon distributions are divided by 3 to 
fit in the plot. 


and growing faster as a; —>• 0. 

The last two points are connected; if there is a deficit in low-x sea quarks of non- 
perturbative origin in our model, this must be compensated by perturbatively generated 
sea quarks from g ^ qq in the PQCD evolution in order to fit HERA F 2 data. This requires 
a large gluon density at small x and a low starting scale Qq, which combined with the mo¬ 
mentum sum rule eq. 0 gives the low width ag of the gluon distribution. If instead of using 
next-to-leading order QCD evolution equations, we use only leading order ones, then the 
evolution is slower which makes the £t to data prefer an even lower starting scale Qq and a 
still more peaked gluon distribution at small x, i.e. a smaller ag, in order to develop a large 
enough qq sea at low x. Thus, these two model parameters depend on the order used in the 
perturbative evolution, whereas the other parameters are not affected much. Furthermore, 
a less good fit to data is obtained with the leading order evolution. 

The Gaussian momentum fluctuations in our model leads to parton number densities 
that go to a constant as x —>■ 0 and hence momentum densities going to zero as x —> 0. The 
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increase of the momentum distributions at x —0 seen in Fig. [TTlis thus caused by PQCD 
evolution from the starting scale Qq, which is forced to be very low. We cannot, however, 
get a sea quark distribution rising as fast as the x“°'^-shape obtained by CTEQ (as seen 
in Fig. ITTbi at their Qq = 1.3 GeV, without an accompanying rise in the xg{x) density. 
This indicates that, besides meson-baryon fluctuations, there is some additional source of 
gg-pairs. One possibility, which seems quite natural, would be a non-perturbative “gluon 
splitting” process similar to the one in PQCD evolution. However, the precise form and 
effect of such a process is not known. Such a contribution wouldn’t have any large effect 
on the asymmetries studied in this work, since the asymmetry data are found at larger x 
and Q^. Therefore we have chosen not to complicate our model by further assumptions and 
parameters to add such a conceivable component. 

IV. CONCLUSIONS 

We have presented a model based on Gaussian momentum fluctuations and hadronic fluc¬ 
tuations of baryons into baryon-meson pairs. The resulting parton momentum distributions 
describes a wealth of data surprisingly well in view of the simplicity of the model. In this 
study, we have focused on different asymmetries of quark distributions in the proton; most 
notably the u-d asymmetry obtained from the forward-backward asymmetry at the 
Tevatron and the d-u asymmetry measured in pp and pd scattering. From the phenomeno¬ 
logical success of the model we can draw some overall conclusions: 

• The Gaussian, based on the law of large numbers, used for the momentum distribution 
may justify a statistical description of the non-perturbative forces conhning partons 
in hadrons. 

• Asymmetries in the quark sea can be explained by hadronic fluctuations, such as 

\p) = ttobo) + |p7r°) -f 

(y-nir |n7r+) 4-... -f- aAi^|AiF+) -4 ... 

From the details of our study, we can make some further assertions: 

• The difference in shape between the u and d valence distributions of the proton corre¬ 
sponds to a difference in the widths of their momentum distributions, which may be 
interpreted as a difference in the spacial region available due to Pauli blocking. 
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• Comparison with the experimental result for the distribution of d{x) — u{x) indicates 
that a larger effective mass is needed for the pion in the |iV7r) states, which may be 
interpreted as an effective way of also including fluctuations with heavier mesons. 

• Besides fluctuations into hadronic states, some other mechanism seems to be needed 
to explain the rise in the sea ^^-distributions as a; —0, used in standard parameteri- 
zations of parton distributions. 

• An asymmetry between the strange and anti-strange nucleon sea is unavoidable, al¬ 
though the magnitude of the asymmetry S~ = dx(xs(x) — xs(x)) cannot be stated 
more precisely than 0.0010 < S~ < 0.0023. This reduces the NuTeV anomaly from 
three to about two standard deviations, leaving no signihcant indication for physics 
beyond the standard model. 

Although our model invoke similar basic hadronic fluctuations as in meson cloud models, 
it differs in important aspects, primarily concerning the basis for the momentum distribu¬ 
tions of both the hadrons in the fluctuations and the partons in the hadrons. As discussed 
above, our model is also more economical in terms of using fewer parameters than even 
the simplest parameterizations of parton densities. Our model provides the parton distri¬ 
butions for all light quarks, as well as the strange sea and gluons, with altogether only four 
parameters giving the x-shapes and three parameters for normalizations. 

Finally, we note that this model could be used to develop a library of parton density 
functions for more general usage. So far, we have not done this since our intention has not 
been high precision functions to be folded with parton level cross-sections in practical cal¬ 
culations. Instead, our primary aim has been to gain understanding of the non-perturbative 
QCD dynamics of the bound state nucleon that is embodied in the parton density functions 
at the low scale Ql, which are otherwise just parameterized. 
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